function [U,U_early,U_late] = fn_J_FB(R1,D0i,lambda,gam,s_min,s_max,dist_s,y1,y2,phi)

c1 = D0i*R1 + y1*D0i;
f  = @(s) util(s.*(fn_rho1(s,phi)-lambda*R1)/(1-lambda)*D0i + y2*D0i,gam).*pdf(dist_s,s);

U_early = util(c1,gam);
U_late  = integral(f,s_min,s_max);

U = lambda*U_early + (1-lambda)*U_late;

end